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Abstract 

Systems composed of distinct complex networks are present in many real-world environments, from so- 
ciety to ecological systems. In the present paper, we propose a network model obtained as a consequence of 
interactions between two species (e.g. predator and prey). Fields are produced and sensed by the individuals, 
defining spatio-temporal patterns which are strongly affected by the attraction intensity between individuals 
from the same species. The dynamical evolution of the system, including the change of individuals between 
different clusters, is investigated by building two complex networks having the individuals as nodes. In the 
first network, the edge weight is given by the Euclidean distance between every two individuals and, in the 
case of the second network, by the amount of time two individuals stay close one another. A third network 
is obtained from the two previous networks whose nodes correspond to the spatially congruent groups. The 
system evolves to an organized state where Gaussian and scale-free-like strength distributions emerge, re- 
spectively, in the predator and prey networks. Such a different connectivity is mainly a consequence of 
preys elimination. Some configurations favor the survival of preys or higher efficiency of predator activity. 

PACS numbers: 87.23.Cc,89.75.Fb,89.75.Hc 
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I. INTRODUCTION 



Despite the impressive development undergone by complex networks along the recent years, 
relatively little attention has been focused on the investigation of systems incorporating more than 
one network, and even lesser attention has been given to the study of specific dynamics (at the 
nodes or of the own network topology) taking place in such networks. Yet, several real- world sys- 
tems are inherently composed by multiple networks. In biological sciences, for instance, protein- 
protein interaction networks are immediately related to transcriptional and metabolic networks yj]. 
Similar entanglements are found in many other areas, such as in the Internet or in the www, which 
are directly related to social and cultural networks, as well as economical constraints [2]. Strictly 
speaking, it is actually hard to think of a real- world system which can be completely represented 
in terms of a single network. 

Some works considering multiple networks have been reported in the literature. The "Solomon 
Network" system was introduced by Erez and collaborators [3]. It consists in a multi-layered 
system with a set of nodes (agents) common to all layers. In each layer, different rules according 
to social and/or economical ties define the interaction between the nodes. On the other hand, the 
coupling between layers is obtained by interaction between variables associated to the respective 
nodes in different layers. Another layered-based formalism was proposed independently by Kurant 
et al. [4]. Using transportation networks, they constructed a logical network representing the traffic 
flows which were mapped onto another network representing the physical infrastructure [|4j, |5|]. 
Park et al. constructed two networks of musicians in which the edges of one network represented 

n 

the similarity while, in the other, they represented the collaboration between the musicians Q6[|. 
The intersection of both networks generated a new dataset where the structural properties showed 
significant differences between both types of networks. We have been particularly interested in 
addressing the dynamics of multiple networks systems. In a previous work [7], we considered 
a system composed of regular and complex networks. While a diffusion pattern evolved in the 
regular network, the complex network were expected to self-organize to control and if possible, 
eliminate the pattern. An interaction rule between both networks was responsible to activate the 
complex network nodes. 

Of special relevance in which concerns multi-layered system are ecological environments 
where the evolution is constrained by the interactions among its components, besides climate con- 
ditions and other external influences. Modelling all such constrains might be difficult and often 
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unnecessary while investigating specific phenomena of interest. Considering an ideal environment 
(favorable climate and abundant supplies), the interaction between animals could be described as 
an interaction function representing the odor intensity [8]. Conversely, the interaction rule could 
be avoided to some extent as in the famous Lotka-Volterra system (predator-prey model) flfloL 
which considers only the amount of individuals of one species as a constrain to the evolution of 
the other species. Adding a diffusive term generalizes the system in which concerns spatial con- 
straints, but does not properly address the question of single interactions Jilt] . Although collective 



behavior may emerge from such simple interaction rules [Q, I3|, LL4J], there are open questions 
related to emergence of spatio-temporal patterns on complex adaptive systems. Since the predator 
prey dynamics involves two types of agents (or individuals) in constant interaction, we considered 
a multi networks view of the system evolution. In this case, each network with moveable nodes 
represents one species, while the interaction rule between the nodes of both networks is defined in 
terms of sensitive fields. 

More specifically, in the present paper we propose a system composed of two interacting 
species, with emphasis on the representation of the spatial position of each individual. One of 
the species (henceforth called predator) relies on the other species (prey) as the only source of 
food. In this scarce food environment, predators move by sensing the presence of preys while the 
preys are expected to sense the predators proximity and move away. Mutual attraction between 
same species individuals tends to generate spatial clusters which imply decrease in the mobility 
of single individuals inside those clusters. Two weighted complex networks are constructed along 
all time steps in terms of the dynamics between predators and preys. One of them expresses the 
Euclidean distance between any two animals as weights. This network is geographical, in the 
sense that the each node incorporate information about the position of the respective individual. In 
the other network, the weights reflect the history of proximity between pairs of animals. In other 
words, the weights in the second network are proportional to the total time each pair of animals 
spend together. A third network is obtained so that each of its nodes corresponds to one spatial 
cluster and the weights provide information about the number of exchanges of animals between 
pairs of respective spatial clusters. 

The current paper begins by reviewing basic concepts related to complex networks and follows 
stating the interaction rules of the proposed model and the parameters of the system. The results 
section discusses the emergence of spatial clusters patterns and then analyse the evolution of struc- 
tural properties of the resulting complex network and their relation with the proposed dynamics. 
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II. COMPLEX NETWORKS CONCEPTS 



A complex network T with TV nodes can be defined by a set V^r) of nodes fa (i — 1,2, . . . , N) 
and a set of edges E(T) connecting a pair of nodes (fa, fa). Given a set W(T) of real val- 
ues, a weight uJ^ hH is assigned to an edge by mapping one element of W(T) to one element 
of E(T). A sub-network k is defined by a set of nodes V(k) C V(r) and a set of edges 
E(k) Q {(fa, fa) : (fa, fa) G £"(r) and fa, fa e V^)}. The sub-network k is connected if and 
only if, any node of V(k) can be reached by any other node of V(k) 



The local connectivity of a node fa can be quantified by its degree fc$, which provides the 
number of nodes directly connected to fa, i.e. the neighbours of fa. When the connections have 
weights, the strength Si of node fa is obtained by summing all weights assigned to connections 
involving fa [18]. The number of closed triangles in the neighbourhood of node fa is given by the 
generalized clustering coefficient cq (eq.Q]) which considers the weights u)^^ of each connection 



established between a node fa and its neighbours fa (20D. 



N,N 
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cc * = Si (ki-1) ^ 2 a ^ a jk (1) 

j,k=l 

where, ay represents the connection between fa and fa such that a^a^a^ = 1 indicates that a 
closed triangle exists in the neighbourhood of node /z^. Finally, the overall structure of a network 
can be summarized by the averaged value of each considered measurement. 



III. INTERACTION MODEL 



A. Movement and Interaction Dynamics 

Consider a squared bi-dimensional region O of size L = 512 in the continuous space Q. We 
randomly distribute TV = 300 animals (or individuals I^, where: i = 1, 2, . . . , N) of each species 
(predator and prey) inside O such that any two species-independent individuals J* and 2j are at 
a minimum Euclidean distance = i? m i n = 4 from each other. At every time step At = 0.1, 
each individual 1 { updates its respective position Pl uVi according to P^. = Pxi, yi + v i^i. The 
resulting sensitive field is understood to provide the velocity vi of each individual at any position 
and time. In such a way, an animal only senses other individuals within its perception radius, 
which is fixed over time as R ner = 128. 
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The interactions between predators and preys are directly related to their species. In case they 
belong to the same species, the sensitive field generated by one over the other is given by a dimen- 
sionless function (eq.[2a|), where the parameters A = 4.0 and /i = 0.04 have respectively units of 
[L] and [L]~ 2 . The alternate characteristic of the function g(rij) is chosen because it provides a 
balance of attraction and repulsion between same species. On the other hand, when the animals 
belong to different species, the individual I { of one species senses Ij of the other species in such 
a way that its reaction intensity is inversely proportional to the proximity between them. This is 
expressed in equation [2bi where the parameter a = 1.0 has units of area [L] 2 . 



hi (2a) 



fj(rij) = -^iij where: ||r^|| < R per (2b) 

The sum over all contributions Ij determines the direction and intensity of an individual J* final 
velocity (eq.[3]). The signs and values of the parameters a and (3 in equation [3] should be selected 
so as to correctly represent the system of interest, and they have units of velocity [£][T] -2 . In 
the present paper, we will investigate different combinations of the parameter j3 for predators and 
preys (eq.0]). 

N-1 N-1 

Predator: a predator <0 and (5 predator >0 

(4) 

Prey: ot prey > and j3 prey > 
Since the model is motivated by an ecological environment, we suppose that preys have plenty 
of food. Consequently, a prey Ij only dies when close enough r^ ?J ^ R e u = 4 to a predator I { . The 
conservation of the number of animals is provided by a feedback process where each eliminated 
prey is promptly replaced by a new one. The feedback process consists in randomly arranging the 
new preys within the region O, while respecting the minimum distance R m in = 4 between two 
individuals as in the initial conditions. As the abundance of preys provides enough resources to 
predators, we suppose that they never die along the simulation interval. 

Although periodic boundary conditions are usually adopted in order to minimize finite size ef- 
fects, an ecological system is often better described by other types of boundary conditions. In fact, 
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animals tend to live inside a specific region O of the ecosystem Q called habitat (Similarly, people 
use to restrict their everyday life within cities or even within home, work and school places). The 
habitat contains not only one species, but also many animal and vegetal species, which consti- 
tute a rich environment. Consequently, the animals are able to interact with a number of species 
enclosed in such a region. The boundaries of the habitat can be related to landmarks such as moun- 
tains and rivers, or even islands. In order to reproduce such real boundary constrains, we adopt 
elastic boundaries conditions (eq. [5]) which restrict the movement region but also allow preys to 
escape while close to the edges. 



B. Complex Networks Construction 

At each time step, a geographical and weighted network fully connected T geo is obtained by 
associating a node to each individual li of both species and setting the weight co^,^ of each 
edge as the Euclidean distance between the respective two individuals. Since the animals are 
constantly moving throughout the sub- space O, the network T geo has a dynamical structure which 
changes at every step. By eliminating the edges above a threshold T = 30 and setting uj^^ = 1 
for the others, we obtain a new complex network Y f geo fragmented in many connected sub-networks 
Kk which we call spatial groups (Fig.QJa). 

The second complex network Y his is constructed by considering the history of proximity be- 
tween every two nodes of T geo (Fig.QJa). The nodes V(T his ) V{Y geo ) are initially fully con- 
nected so that each edge is set with = 0. Whenever two individuals 1 { and Ij fall close 
enough (nj ^ T = 30), the weight w 1 *^ = + 1 is updated. The connections involving 
dead preys are immediately eliminated. 

The resulting network T res is obtained by a merging mechanism which is done by associating 
a node v k to each sub-network K k of Y' geo such that a new set of nodes V(Y res ) is obtained. If 
M is the number of edges between the sub-network K k and K\, the weight E^i (eq. [6]) of the 
corresponding edge (v kl vf) in T res is given by the average value of established between the 
members and jij of different spatial groups K k (Fig.QJb). 



If x(y) > L then x{y) <— 2L — x{y) 
If x(y) < then x(y) <— \x(y)\ 



(5) 
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Figure 1: (a) A sample of both networks where the nodes represent the individuals. Light-gray edges belong 
to T f geo and black edges to network r^ s . Dashed lines emphasize the sub-networks k^. (b) The resulting 
network Y res is obtained by merging V geo and r^ s . Each node of r res is related to one sub-network Kk 
in Figure [U-a. In both Figures, the edge thickness is proportional to the respective weight. 



IV. CLUSTER EMERGENCE 
A. Spatio-temporal patterns 

The interaction between different species, along with the prey replacement mechanism, are fun- 
damental for providing a non- stationary state. Spatio-temporal cluster patterns uM emerge from 
the dynamics due to interaction between same specie individuals. In the absence of one species 
the system evolves to a stable state with many small static clusters (Fig. [2]). Contrariwise, in the 
presence of interacting preys and predators, but null attraction between same specie individuals, 
the system evolves to a non-stable state, without emergence of clusters (Fig. [3]). 

In order to take into account the diversified behaviours observed in nature, we now investigate 
the emergence of patterns considering different values of j3 for each species (/? = 0.5, 0.125 and 
0.0625). We start by considering null interaction between preys while the predators are able to feel 
their counterparts ((3 pr edator > and f5 prey = 0). Observing the evolution of the system when three 




(6) 
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(a) t = 250000 (b) t = 500000 
Figure 2: Snapshots of the emergent spatio-temporal pattern for one single species with j3 — 0.125. 



(a) t = 250000 (b) t = 500000 

Figure 3: Snapshots of predators (red - on-line) and preys (black - on-line) spatio-temporal patterns with 
null attraction (f3 = 0) between same species individuals. 

distinct values of (3 pr edator and a = 32 are considered (Fig. H]), we clearly identify the emergence 
of predator clusters which differ according to the f3 predator intensity. 

When the attraction between predators is stronger (i.e. f3 predator = 0.5, see Figure @]-i), non- 
uniform clusters emerge in a few time steps (Fig.|4]-i,a). Mutual attraction between different clus- 
ters progressively generates larger clusters, with the smaller clusters tending to be attracted by the 
larger ones (e.g., as is the case with the two clusters in the middle of Figure @]-i,a). Although the 
clusters present several shapes in the first steps, they evolve towards spherical shapes because of 
symmetrical internal forces and increase in the number of members (Fig.@]-i,b). After a transient, 
the number of clusters seems to stabilize. The smallest clusters can move throughout the space 
following the concentration of preys (Fig.|4]-i,d) while the denser clusters repel strongly the preys. 
Consequently, the cluster net velocity becomes slower than that of the preys, allowing the latter 
to escape. As a consequence of the boundary conditions, such dense clusters tend to concentrate 
in the central region of O while the preys move toward the periphery (Fig. @]-i,e). Therefore, we 
expect a resulting giant cluster of predators to appear after a long period of time. 
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50000 


(b) t = 100000 


(c)t 


= 200000 
(iii) 


(d) t = 300000 


(e)t = 


500000 



(a) t = 100000 (b) t = 200000 (c) t = 300000 (d) t = 400000 (e) t = 500000 

Figure 4: Snapshots showing the pattern evolution of predators (red - on-line) and preys (black - on-line) 
when the attraction between predators is given by (i) (3 pre dator = 0.5, (ii) (3 pre dator = 0.125 and (iii) 

/^predator = 0.0625. 

Emergence of clusters also occur in the second configuration with (3 pr edator = 0.125 and 
ftprey = (Fig. @]-ii). However, they are visually more uniformly distributed (Fig. @]-ii,a) when 
compared with the previous configuration (Fig. @]-i,a). Since the attraction intensity is smaller, 
dense clusters emerge later along the dynamics (Fig. @]-ii,c) while the coarser clusters tend to be 
maintained along time (Fig.[4]-ii,d). The last obtained time step (Fig. Sl-ii,e) resembles the initial 
stage of the first configuration (Fig. @]-i,a), suggesting that the system would possibly evolve to a 
similar state. However, the giant cluster is not expected to emerge because of the weaker attraction 
between predators in this case. 
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(a) t = 10000 (b) t = 75000 (c) t = 150000 (d) t = 500000 




(a) t = 50000 (b) t = 100000 (c) t = 200000 (d) t = 500000 

Figure 5: Snapshots showing the evolving spatio-temporal patterns of predators (red - on-line) and preys 
(black - on-line). The attraction between same species individuals is given by (i) (3 = 0.5 and (ii) (3 = 
0.125. 



The last snapshots sequence (Fig. @]-iii) presents the evolution when the lowest attraction in- 
tensity between predators is considered ((3 pre d a tor = 0.0625 and (3 prey = 0). In this configuration, 
we cannot identify dense clusters up to the last time step. However, spatio-temporal patterns of 
uniformly distributed small clusters emerge and evolve in a non- stationary way with predators 
and preys moving close to each other. The pattern consists of clusters with different sizes and 
shapes, indicating that predators are able to move throughout sub-space O, moving between clus- 
ters (Figs. @]-iii,d and@]-iii,e). Any attempt of escape by preys is promptly checked by predators, 
which can easily enclose any group of preys since the predators are not strongly attracted in this 
configuration. 

When we consider null attraction between predators and different attraction intensities between 
preys, the spatio-temporal pattern tends to include uniformly distributed predators and preys (sim- 
ilar to Figure [3]). Actually, small clusters of preys emerge but are readily eliminated because of 
the strong predators attraction, suggesting a cyclic process. Although we can identify some larger 
clusters with higher attraction intensity, the pattern is nearly independent of (3 prey . 
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The last experiment consists in introducing attraction between individuals of both species, i.e., 
both predators and preys are able to feel their counterparts ((3 predator = f3 prey — (3 > 0) while they 
can also feel each other (i.e. a = 32). The high level of attraction j3 = 0.5 resulted a pattern 
of dense clusters of predators while smaller clusters of preys have also been observed along the 
first time steps (Fig. [5]-i,a). The simultaneous attraction between members of each species gener- 
ated distinct regions occupied by different species (Figs.[5]-i,b to[5]-i,d). The higher concentration 
of preys in some clusters increased the intensity of the sensitive field such that the clusters of 
predators could move around a larger region (following the prey clusters) when compared with 
the configuration with null attraction between preys (Fig. [4]). The clusters of preys often disappear 
and emerge as a result of this attraction (Fig. (5]-i). On the other hand, the clusters of predators 
become denser over time because of the non-uniformity in the distribution of the clusters of preys 
(Fig.[5]-i,d). Since the preys are non-uniformly distributed in space, they create attractor regions. 

The snapshots from Figure EJ-ii show the evolution of the system with lower attraction 
(3 = 0.125 between same species individuals. In this case, the clusters are not so dense and they 
are more uniformly spatially distributed (Fig. [5]-ii,a). The preys are also organized in clusters of 
smaller sizes (Fig.[5]-ii,b). The sequence of snapshots suggests a dynamical behaviour where some 
clusters aggregate more predators while the greatest part of clusters divide and recombine in new 
sets of animals (Fig. [5]-ii,c). The clusters of predators follow the concentration of preys, which 
become higher in some regions over time. 

In the last configuration where the attraction has the smallest value (/3 = 0.0625), the density 
of clusters in the pattern was smaller than in the other cases (Figs.EJ-i and[5]-ii). Actually, the situ- 
ation is similar to the absence of attraction between preys and small attraction between predators 
(Fig.|4]-iii). However, in this case we could identify some small clusters of preys which constantly 
appeared and disappeared along time, while the clusters or predators resulted very small. 

B. Prey Elimination 

Spatial cluster structure emerges in animal behaviour as a consequence of several factors. Usu- 
ally, group structure is a consequence of some type of similarity which provides protection against 
external agents. Protection not only against other individuals, but also against ideas and other 
cultural behaviour. In our model, the cluster structure depends on the attraction intensity between 
same species individuals and can favor species according to the attraction configuration. Although 
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(a) ^predator > and f3 prey = (b) ^predator = and f3 prey > (c) ^predator > and f3 pre y > 



Figure 6: Evolution of the amount of eliminated preys (average and standard deviations). We considered 
three intensities of attraction between same species animals. 



the preys death by itself is not necessarily the final purpose of our model, the number of elimi- 
nated preys indicates the importance of clusters and how they constrain the spatial movement of 
predators and preys. For instance, attraction between preys increases the sensitive field of preda- 
tors because of the larger number of preys close one another. Consequently, the predators velocity 
directed to the prey clusters increases while the cluster structure constrains the movement of single 
preys. Contrariwise, the attraction between predators generates strong fields which allow preys to 
escape faster than the predators can move. This situation results in lower rates of preys deaths 
(Fig.[6]-a), while preys were quickly eliminated in the previous case (Fig.[6]-b). 

The preys death rate is linear in almost all configurations, except in the cases with /3 pr edator — 
0.5, where the growth rate is faster along the initial 200000 time steps, slowing afterwards (Figs. [6}- 
a and[6]-c). This two-slope behaviour is a direct consequence of the emergence of dense predator 
clusters about this time step, decreasing the preys death rate. The existence of attraction in both 
species (Fig.[6]-c) has similar effect as when attraction is allowed only between predators (Fig.[6]-a). 
However, the cumulative number of eliminated preys is larger in the latter case (Fig. [6]-a). These 
results suggest that the spatial organization of animals has crucial importance for their survival. 
Therefore, according to the predator species, there is no difference for preys to be organized in 
clusters or not. Since the hunting ability of predators is similar, the preys organization should be 
determined by other environmental reasons. In this sense, the choice of preys by predators might 
be directly related to the way in which the preys are spatially organized into groups. 
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V. STRUCTURAL PROPERTIES 



The mechanism proposed to build the complex networks generates bi-partite networks with 
two types of nodes, respectively: predators and preys. Eliminating the connections between nodes 
of different types, two sub-networks are obtained which are henceforth named predator Y predator 
and prey T prey networks. We shall analyse the enrollment and movement of individuals between 
spatial clusters by investigating local structural properties of both resulting networks. 

A. Evolution of the structural properties 

Considering the predator network T predator . In the absence of attraction between preys 
(ftprey = 0), the average degree (k) increases quickly within the 150000 time steps, independently 
of the attraction intensity between predators (Fig.[7]-i). After this stage, the average degree nearly 
stabilizes about (k) = 8.5 with f3 predator = 0.5 (Fig.[7]-i,a). However, it exhibits a slower increase 
in the other configurations (Figs. [7]-i,b and[7]-i,c), suggesting that stabilization will also be even- 
tually reached. Actually, because of the movements of predators between the spatial clusters, this 
measurement should stabilize after the network becomes completely connected. This will happen 
at (k) ~ 40 in cases of j3 predator = 0.125 (Fig. [7]-i,b) and of f3 predator = 0.0625 (Fig.[7]-i,c), which 
corresponds nearly to the average number of spatial groups or number of nodes in the respective 
network T predator (the average value is obtained after an initial transient). Because the number of 
groups is smaller when l3 predator is larger, the threshold becomes smaller in that case (Fig. [7]-i,a) 
and it possibly decrease whether the system converges to a giant component of predators. 

In the proposed growth method, the node strength relates to the history of proximity between a 
node and the members of its group, and inversely, to the rate of members exchanges between two 
spatial groups. The average strength (s) evolution seems to take place in stages when the attraction 
intensity between predators is larger (f3 predator = 0.5 - Fig.[7]-ii,a). Possibly, the stages in which the 
growth rate is slower correspond to the intervals where the exchange rate between groups is small 
(Fig.[7]-ii,a). Since the predators spend much time together, the weight of connection (/i^, jij) 
in the network T his increases considerably and when two or more groups merge, the average 
strength suddenly rises. A lower attraction intensity implies on more circulation of predators in 
space and, consequently, constant association and disassociation to several groups. The level of 
proximity between predators tends to be smaller in the configuration with (3 predator = 0.0625 
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(i) 




Time step Time step Time step 

(a) f3p re dator = 0.5 (b) flpredator = 0.125 (c) ^predator = 0.0625 

Figure 7: Evolution of the average structural properties in T pre d a tor when three intensities of attraction 
between predators and null attraction between preys are considered, (i) Average degree (k) Average 
strength (s) and (iii) Average clustering coefficient (cc). The standard deviation is shown to 50% of the 
original values in cases (i) and (ii). 

(Fig. |7]-ii,c) when compared to the other two cases (Figs. |7]-ii,a and |7]-ii,b), where the predators 
movement is more constrained. As a consequence, in the last two cases, the predators spend 
more time together with the same partners. The larger strength in case of (5p redator = 0.125 is a 
consequence of the larger number of spatial groups (nodes of Y res ). 

A smaller attraction intensity between predators implies in a faster increase in the average 
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(c) 



Figure 8: Evolution of T pre d a tor average structural properties when the attraction between same specie 
individuals is given by (3 = 0.5. The standard deviation is shown to 50% of the original values in (a) and 
(b). 



clustering coefficient (cc) within the first time steps (Fig.[7]-iii). In the second growing stage (from 
nearly 25000 to 250000 time steps), the growth rate is higher for f3 predator = 0.5 (Fig. [7]-iii,a). 
The latter effect is possibly a result of the decrease in the number of spatial groups in the first 
configuration (Fig.|4]-i). Since the clustering coefficient measures the local connectivity (between 
common neighbours of a reference node), and the growth of (cc) (Fig.d-iii) is higher than the 
growth of the corresponding (k) (Fig. |7]-i), the connections occur mainly between nearby spatial 
groups. In other words, the movements required for following preys is considerable slow. A 
predator moves within a small Euclidean distance and then associates to one spatial group; after 
a while, the same or other predator moves again to another close spatial group in a non-stationary 
way. 

Considering attraction between preys {(3 prey > 0) and null attraction between predators 
(Ppredator = 0), the structural properties of the evolution is similar to the situation with [3 predator = 
0.0625 and f3 prey = 0, except for the configuration with f3 prey = 0.5. In that case, the stronger 
attraction of predators by dense clusters of preys implies on more movements of predators and, 
therefore, more connectivity between the nodes in comparison to the other two configurations. 
The absence of attraction between predators and the higher concentration of preys due to mutual 
attraction produces an overall behaviour where the predators have freedom to move throughout 
sub-space O. As a consequence, the average measurements values are slightly higher in this case 
in comparison to the situation with small attraction between predators (similarly to configuration 
of Figure |7]-c), although presenting the same stages of evolution. 
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In the last experiment, involving attraction between predators {j3 predator > 0) and between preys 
(Pprey > 0), the evolution of the average degree shows again an apparent stabilization in all cases. 
The number of nodes in the resulting network T predator (or number of connected components k) 
is about TV ~ 9 in the last time step observed (after a large transient). Since TV is higher than 
the plateau in the graph of Figure [8]-a, the resulting network is not completely connected yet. 
We verified that at the maximum value (between 150000 and 200000 time steps), the number 
of spatial groups is approximately 1.5 of the average degree value. Although the node average 
connectivity has decreased after this interval, the value increased if compared with the number 
of nodes in the network. This effect is related to the merging of clusters, which happens more 
frequently in the configuration where the attraction is higher (Fig.[5]-i). This effect is not observed 
in the other two configurations, where the average number of groups is nearly constant (about 
45) over time and the curves are similar to the case in Figure [7]-i. Comparing the number of 
nodes and the stabilization threshold, we observed that the resulting network is almost completely 
connected in the last time steps. The average connectivity among the nodes grows faster in the 
configuration with the smallest value of (3 as a consequence of higher mobility. The average 
strength presented evolution similar to the first experiment (Fig.[7]-ii). The main difference regards 
the stronger attraction between preys which generated clusters of preys. The non-existence of 
uniformly distributed preys implied in less movements between predator clusters due to absence 
of stimuli. Consequently, the strength values are higher (Fig. [8]-b) when compared with the first 
experiment. 

The decrease of the mobility between spatial groups is also observed by comparing the evo- 
lution of the clustering coefficient in Figure [U-c with Figure [7]-iii,a. Once again, the clustering 
coefficient growth rate is very high within the first time steps (up to 10000), maintaining a nearly 
constant growth henceforth. Actually, the evolution seems to stabilize after t = 450000 in all 
configurations. Differently from the first experiment, the system may not evolve to a state where 
all nodes become connected. The effect is more probable in the configuration with (5 — 0.5 pro- 
vided the clusters do not merge in a single giant component. The movement of nodes between 
two clusters become less probable after a certain time step. This is a consequence of the fact that 
each cluster now has higher density and is considerably far from other clusters. The absence of 
stimuli because of prey clusters emergence and the mutual attraction inside the predator clusters 
also contributes to such a decrease of movements. 
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B. Strength Distribution 



The average measurements provide global information about the complex network. However, 
the network internal structure can be completely different even when presenting similar average 
values J 1811. The internal structure is responsible to constrain dynamical processes (e.g., cascade 



failures 



210 or epidemics 



Q, 



220) in the network. Consequently, it is fundamental to 



characterize the resulting network in terms of measurements related to their internal topological 
properties. 

The strength distributions (histograms) have different shapes in the three resulting networks 
(Fig. [9]) when we consider attraction between predators ((3 pre dator > 0) and null attraction between 
preys (f3 prey = 0). Although the configuration with (3 predator = 0.5 (Fig.[9]-a) presents undefined 
shape, we see that the maximum strength increases considerable between 100000 and 300000 time 
steps. The number of nodes with large strength values also increases at the next time steps consid- 
ered (Fig.[9]-iii,a), indicating that the mobility between groups decreases with time, which can be 
seen by the existence of nearly stationary (in shape and number of members), dense spatial clusters 
in FigureS Contrariwise, the other two configurations (Figs.[9]-b and[9]-c) present a Gaussian-like 
strength distribution with a characteristic scale given by the average value (the average can be ver- 
ified by comparison with graphs of Figure [7]-ii). Although presenting a characteristic membership 
permanence time, some individuals constantly change between their respective groups (decreasing 
the strength), while others associate in a group for a long time (increasing the strength) and barely 
move between different groups. The predators in the last case can be seen as the core of the group, 
since they spent much time together and are responsible for maintaining the group united. 

In the case of the prey network T prey , the strength distribution (Fig. [TOj) resulted in a com- 
pletely different function when compared to the predator network (Fig. (9]-iii). In this case the 
histograms remind a power law function (especially along the largest values) with the maximum 
strength value being larger in case of larger l3 predator (Fig. QUa) and the other two configurations 
apparently presenting two-slopes (Figs.[K)]-b and[K)]-c). The shape of the function suggests that 
the system converged to an organized state. In the last time step considered, few preys stayed to- 
gether for a long time (higher strength), while their majority moved away or died frequently (lower 
strength), such that the emergent distribution is scale-free-like (Observe the high number of null 
strength nodes). Although a low strength value could indicate a high rate of changes between 
groups, the low values observed in the preys network is mainly because of preys deaths since the 
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Figure 9: Strength distribution for predator networks at three time steps: (i) t = 100000, (ii) t = 300000 
and (iii) t = 500000. We consider three attraction intensities between same-specie individuals. 



maximum strength value is much larger in the predators network and the amount of eliminated 
preys decreases with larger f3 predator (Fig. [6|-a), while the maximum strength value increases for 
larger (3 predator (Fig.[H). 

Null attraction between predators (P pre dator = 0) and attraction between preys ((3 prey > 0) 
generated a similar pattern as in the case with l3 predator = 0.125 and (5 prey = 0, though exhibiting 
a Gaussian-like strength distribution. Comparing with the case when only attraction between 
predators is considered (Fig. [9]-c), we verified that the maximum strength is considerably higher 
at the same time step observed in that case. It suggests more movements between groups in the 
configuration with (3 prey > and (5 predator = because of the stronger attraction provided by the 
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Figure 10: Strength distribution for T prey at t = 500000 when three intensities of attraction are considered 
between predators. 



high concentration of preys. Once again, the prey elimination mechanism implies in short life-time 
for preys and consequently, suggesting a power law distribution along the highest values. 

In the last experiment, considering attraction between same species individuals, the strength 
distributions present some interesting effects. The configuration with larger j3 converged to a 
power law-like distribution (Fig. Qjja) due to formation of dense prey clusters which implied in 
faster predator clusters emergence in some regions of the space. Dense predator clusters have 
strong attraction fields and therefore, any predator close to them is readily attracted and barely 
escapes from the cluster. Consequently, only small groups change individuals frequently. In case 
of j3 = 0.125, the system evolves from a Poisson-like state to a combination of Gaussian-like on 
the left with a power-law-like on the right (Fig.QT]-b). At last, (5 = 0.0625 (Fig.[]T]-c) is similar 
to the cases with small (Fig. [9]-c) or null attraction between predators, suggesting that a reason- 
ably weak attraction intensity will result in the same non- stationary state. In other words, weak 
attraction between predators and/or between preys does not produce different collective behaviour 
in the system. 

The power law-like distribution is observed in the prey network T prey for large values as in the 
other configurations (Fig.O. The effect confirmed the hypothesis that such a scale-free structure 
emerged because of the elimination of preys and not because of the attraction mechanism be- 
tween two individuals. The longer life-time of preys is identified by the larger maximum strength 
observed in this case in comparison to the other configurations (Tig.fTOl). 
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Table I: Summary of the main properties (spatial and topological) according to different attraction intensities 
between same species individuals. 
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Figure 11: Strength distribution for T pre( i ator at t = 500000 when three attraction intensities between same 
specie individuals are considered. 
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Figure 12: Strength distribution for T prey at t = 500000 with three attraction intensities between same 
species individuals. 



C. Clustering Coefficient Distribution 

The clustering coefficients have nearly the same distribution in all considered configurations 
for the predator network (Fig. [13]). The only difference is observed at the scale of the distributions 
since smaller (3 pre dator implies more movement and, consequently, faster connections between 
groups. The normalized distributions become narrower and present characteristic scales which 
move to the right over time, indicating that all groups will become connected after a long period 
of time. The shapes of the curves indicate that higher attraction between predators constrains 
the movement and postpones creation of triangles between common neighbours of a reference 
node. In case of weaker attraction, predators are able to move between groups and consequently, 
the amount of nodes with larger clustering coefficient becomes higher (Fig. [T3]-c). Although the 
number of nodes with higher clustering coefficient constantly increases over time, there are many 
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Figure 13: Clustering coefficient distribution for T pre( i ator at three time steps: (i) t = 100000, (ii) t 
300000 and (iii) t = 500000. We consider three attraction intensities between predators. 



nodes with small values of cc for larger (3 predator (Tig.[T3la). 

In the preys network T prey , all configurations resulted in a clustering coefficient distribution 
with a clear division (Fig. [14]), which suggests the existence of two distinct states in the system 
(nearly constant over time). A large number of preys has small connectivity in their neighbourhood 
while a small number presents a high level of local connectivity among its neighbours. A small 
clustering coefficient is a consequence of two effects: complete absence of movement or a high 
death rate, which avoid triangle formation. Since the frequency of cc null values increases as we 
decrease the attraction between predators (Fig. [14]), we conclude that the main fact behind this 
effect is exactly the death of preys. As shown before, such deaths increase as Ppredator decreases 
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Figure 14: Clustering coefficient distribution at t = 500000 for T prey when three intensities of attraction 
between predators are considered. 



(Fig. [6]-a). In case of attraction between preys, the effect is the opposite, i.e. the number of dead 
preys and of nodes (predator groups) with null cc increases with f3 prey . Since the configuration with 
fiprey = 0.5 has regions with high concentration of preys, which are not completely static and are 
not eliminated frequently too, the probability to establish closed triangles increases and explains 
the higher frequency observed at the maximum value of cc (cc = 1). The existence of groups 
with maximum clustering coefficient indicates that some preys stay a long time in the system 
and, although they are eliminated, other preys keep some spatial groups united and permit triangle 
formation (about 10% of the nodes in T prey ). With attraction between same species individuals, we 
observed a higher frequency of intermediate clustering coefficient values. The existence of clusters 
in both species resulted in extended life-time for the preys, which allowed more movement and 
connections between nearby spatial groups. 



VI. CONCLUSIONS 



The current paper has focused on multiple networks systems involving interacting species. 
Since the system was originally motivated by ecology, we assumed that preys are eliminated when 
they move sufficiently close to a predator, but new preys are randomly displaced in the system in 
order to replace those which are consumed. Spatio-temporal clusters emerge from the dynamics, 
which involves movements of individual between clusters. Two complex networks containing the 
same set of nodes are considered. In one network, the weights represented the Euclidean distances 
between two individuals. In the other network, the number of steps during which two individuals 
were close enough are considered as weights. By merging both networks, we obtained a third 
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complex network whose nodes represented spatial groups defined by the connected sub-networks. 
Such a growth mechanism implied the connections to incorporate information about the history of 
the individuals movement. 

Several configurations, defined in terms of the intensity of the attractions between preys and 
predators, were considered in our simulations. The increase in the attraction between same species 
individuals generated dense clusters whose sizes and shapes were a result of the same species 
attraction intensity as well as the other species attraction features. Dense predator clusters con- 
strained the movement of single predators, while dense prey clusters generated strong fields and 
attracted many predators. As a consequence, the rate of eliminated preys was higher when the 
preys were organized into clusters and predators had null attraction between them. We have ob- 
served that group organization could benefit the species according to the intensity of intra- species 
attraction. On the other hand, some configurations (with small attraction) evolved to states where 
the group organization had few or no advantage to any of the two species. 

By using the complex network theory, we observed that the average degree of the resulting 
predator network increased up to a threshold corresponding to the number of spatial groups, i.e., 
the system evolved to a fully connected state. Since the average clustering coefficient increased 
faster than the average degree, we concluded that members exchanges between nearby spatial 
groups were more frequent than between far away groups. The analysis at the last time step showed 
that the system converged to an organized state. The predator network presented a Gaussian-like 
strength and clustering coefficient distributions in nearly all configurations while approximately 
scale-free strength and a polarized clustering coefficient distributions emerged in the case of the 
preys network. The prey elimination mechanism was responsible for generating the observed 
structure in the prey network. Since preys are eliminated frequently, they are not able to move 
long enough throughout space and establish connections between all clusters, as was the case with 
the predators. 

In order to further investigate the model, we propose the following future developments: (i) 
study of the effect of individuals density and other scale properties of the system, (ii) investigation 
of clusters emergence when considering other attraction functions and more species interacting 
together, (iii) analysis of collective phenomena (e.g., epidemics and opinion formation as a result 
of dynamical propagation of diseases and ideas from specific individuals), and finally (iv) appli- 
cation of the same methodology in empirical data, for example, to study the movement of animals 
in the field. 
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